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Nuclear magnetic resonance in the heavy fermion 
superconductors 



1. Introduction 

Nuclear Magnetic Resonance (NMR) has impacted several branches of science, not the 
least of which is the study of strongly correlated electrons in condensed matter. The 
pioneering works of Rebel and Shchter [H [2] and Masuda and Redfield |3] revealed the 
power of this technique to the condensed matter community. These NMR experiments 
measuring the spin lattice relaxation rate in the superconducting state of Al showed 
the presence of a superconducting gap, and more importantly that the Cooper pairs are 
a coherent superposition of two particles. Their results clearly demonstrated that the 
two-fluid picture of superconductivity, consisting of single-particle states, was untenable. 
Indeed, in their seminal paper on the theory of superconductivity, Bardeen Cooper and 
Schrieffer point out the role of both NMR experiments and ultrasonic attenuation in 
support of their theory [1]. 

NMR played a central role in the study of the high-Tc cuprates, and continues to 
yield important information about novel superconductors [3 [6] . In the last decade or so, 
there has been a renaissance in the study of the heavy fermion superconductors. These 
materials were first discovered in 1979, but were largely overshadowed by the tremendous 
effort focused on high-Tc materials [Tl[8],[9]. Recent discoveries of the CenMmIn3n+2m 
family of compounds, various actinide-based superconductors, the breakdown of Fermi- 
liquid theory, and the emergence of superconductivity in the vicinity of quantum 
phase transitions have turned the attention of the condensed matter community to 
the nature of the superconductivity, magnetism, and new broken symmetries that 
emerge in these compounds. This article is intended to briefly review some of the 
important contributions of NMR to the study of these phenomena, with an emphasis on 
a discussion of NMR as a general technique for investigating superconductivity rather 
than a comprehensive review of the literature. In particular, this article is directed at 
non-NMR specialists who wish to understand and critically evaluate new data. For an 
in-depth review of NMR in heavy fermions, the reader is referred to Ref. |10j . 

Upon reflection, one might expect a priori that the Meissner effect in a 
superconductor should preclude any NMR experiments. Indeed, many NMR 
experiments in superconductors are dominated by the response of the superconductor 
to magnetic fields. NMR studies of superconductors can be divided into two categories: 
those which investigate processes driven by (i) hyperfine fields and by (ii) orbital 
supercurrents. Hyperfine fields arise from interactions between the nuclear and electron 
spins. In metals and superconductors, the dominant contribution to relaxation of 
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the nuclear spins usually arises from scattering from the quasiparticles. Orbital 
(super) currents arise in superconductors as a response to screen external magnetic 
fields. Typically these currents are manifest in a vortex lattice, and give rise to an 
inhomogeneous magnetic induction. Since the nuclear spins couple to this magnetic 
induction through the Zeeman interaction, the NMR spectrum is a sensitive probe of 
the physics of the vortex lattice. 

In order to fully understand the behavior of nuclei in a superconductor, it is 
crucial to first understand their behavior in conventional metals. The principles of 
NMR in metallic systems is discussed in sections ^ and ^ Section ^ addresses the 
hyperfine interaction in superconductors, which probes the low energy properties of 
the superconducting gap function. Section ^ focuses on measurements that couple 
to the orbital supercurrents, and reveals physics of the vortex lattice. In section ^ we 
review NMR measurements in the superconducting state of the heavy fermion materials. 
NMR measurements in the high temperature superconductors has been instrumental to 
elucidate the physics not only of the superconducting state, but the low energy physics of 
the spin system. Aside for a brief discussion in §l]on the measurement of the symmetry of 
the superconducting order parameter in the cuprates via spin echo decay measurements, 
this article will not attempt to review the vast literature of NMR in the high temperature 
superconductors. We refer the reader to [5] for details. 

2. NMR Basics 

2.1. Hamiltonians and Spectra 

Nuclear spins in condensed matter constitute are a nearly ideal example of an isolated 
ensemble of coupled particles that experience a weak coupling to an external bath. This 
weak coupling, the nuclear-electron couphng, gives rise to small perturbations of the 
nuclear spin Hamiltonian, as well as provides a channel for relaxation and decoherence of 
the nuclear spin system. It is via this interaction that the nuclei can provide information 
about the electronic degrees of freedom in strongly correlated electron systems. The 
nuclear spin Hamiltonian in correlated systems can be written as the sum of four 
interactions: 

H = Hz + + T^dip + ^e-n- (1) 

The first term describes the Zeeman interaction between the nuclear spin I and the local 
magnetic field, Hq: 

Hz = ih^o ■ i, (2) 
where the gyromagnetic ratio 7 is related to the nuclear moment /i by 7^ = 
g^J'N^/ I{I + !))• Here /itv is the nuclear magneton, and g is the nuclear g-factor. This 
term is usually the dominant interaction for nuclei in an external field, with 7^/9 ~ 100 
MHz. Since h/ks ^ 48yuK/MHz, the nuclear spin system can be treated in the high 
temperature limit, as most phenomena of interest in correlated electrons are manifest 
at temperatures > 1 mK. 
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The second term in Eq. ([T]), the quadrupolar interaction or the charge part of the 
electron-nuclear interaction, is given by: 

= iT^^f^P^' - + Vill - I'y)], (3) 

where e is the electron charge, Q is the quadrupolar moment, V^^ are the components 
of the EFG tensor, and rj = {Vxx — Vyy)/Vzz is the asymmetry parameter of the 
EFG. The NQR frequency (in zero field) is defined as: uq = z/^-^/l + r/^/S, where 
i^z = eQ/2I{2I — l)h. For a site with axial symmetry [rj = 0), there are 2/ — 1 
quadrupolar resonances at frequencies nuq, where n = 1, . . . , 2/ — 1. If rj > 0, then the 
resonances are not equally spaced. The EFG is fully characterized by three parameters: 
z/2,?7 and a unit vector, q, in the direction of the principle axis of the EFG with the 
largest eigenvalue. The quadrupolar interaction vanishes for / = 1/2, but can be 
significant for nuclei with sufficiently large quadrupolar moments and / > 1/2. Typically 
{T^q) /h ~ 1 — 100 MHz. In practice, this means that the degeneracy of the nuclear spin 
manifold is lifted even in zero field, and one can detect resonances in the absence of 
an applied field. This technique is known as Nuclear Quadrupolar Resonance (NQR), 
whereas the term NMR is reserved for resonance experiments in field. This is particularly 
important for superconductors, since it allows one to measure the spin lattice relaxation 
in zero field, without the infiuence of the vortex lattice. 

The third interaction in Eq. ([1]) is the dipolar interaction between nuclear spins in 
a lattice, and is given by: 

^ . ^ 1 [ N ■ /^fc _ 3(/ij ■ rjfc)(/ifc ■ Yjk) \ 

where Yjk is the distance between the nuclei moments j and k with moments [i = jhl. 
In most solids, the dipolar interaction gives rise to a static broadening of the lineshape, 
and in most cases the lineshape can be well approximated by a Gaussian with second 
moment given by: 

j,k i'^ 

Usually (Aa;^) is temperature independent, and depends on the structural details of the 
particular compound of interest. Typically (Acu^) ~ 1 — 10 G^. For studies of vortex 
physics in superconductors, it is preferable that (Acu^) is as small as possible, since the 
inhomogeneous field distribution in the mixed state of a type II superconductor broadens 
the lineshape in a quantifiable fashion. Since this vortex lattice lineshape broadening 
is typically only on the order of a few Gauss, one is often faced with the fact that the 
intrinsic dipolar broadening limits the spectral resolution (see ^1). Furthermore, the 
quadrupolar interaction ([3]) can also contribute to line broadening, particularly when 
there are lattice strains which can give rise to distributions of EFGs. 
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The fourth term in Eq. ([T]) is the magnetic part of the nuclear-electron interaction, 
and is given by [llj: 



where S, and Ij are the spin and angular momentum (with respect to the nucleus) of 
the i^^ electron. The first term of Eq. ([6]) is the orbital interaction, which gives rise to 
the chemical shift in materials where orbital angular momentum is unquenched, as well 
as in superconductors with supercurrent distributions. The remaining three terms are 
the dipolar and contact part of the hyperfine interaction, which is usually the dominant 
relaxation mechanism in metals and superconductors. Both terms can be understood 
in terms of extra fields experienced by the nucleus, either by a hyperfine field, Hhf, 
or by a magnetization NL^rb created by a current distribution. The former typically 
is the primary interaction between nuclei and quasiparticles in superconductors, and 
the latter couples the nuclei to the vortex lattice in superconductors. In most strongly 
correlated electron systems, the hyperfine interaction is complicated by the presence of 
transferred hyperfine interactions, given by a complex exchange process involving orbital 
overlaps between different ions in a material. In practice, the hyperfine interaction is 
approximated as: 



where A is an on-site hyperfine coupling (typically dominated by core-polarization 
effects), and Bi are transferred hyperfine coupling to the neighboring spins, S(rj). 
These couplings are usually taken as empirical values, since accurate calculations of 
these values involves detailed electronic structure calculations, which are not always 
available, especially in strongly correlated systems [12]. Eq. ([7]) is the basis of the Mila- 
Rice-Shastry Hamiltonian in the cuprates [HI |14] , and is a good approximation to the 
interactions in most heavy electron systems |15j . 

2.2. CW versus pulse techniques 

Equation ([1]) lifts the degeneracy between the nuclear levels, and the spectral 
properties can be understood in terms of the energy level differences between the 
eigenvalues of Ti. A basic NMR spectrometer measures the energy absorbed by an 
ensemble of nuclear spins as a function of either fixed field while frequency is swept, 
or vice versa. Such an experiment describes a basic continuous-wave (CW) spectrum 
measurement. Over the last three decades pulsed techniques have been developed such 
that a nuclear spin system can be excited by a single short pulse (of duration ^1 — 10 
/isec). The spectrum is then given by the Fourier transform of the response to the 
pulsed excitation. A spin echo is the response of a nuclear spin system to a series of two 
pulses. The Fourier transform of a spin echo is also identical to the spectrum, with the 
advantage that the response of the system is separated temporally from the excitation 
pulse, eliminating a source of noise. 






(7) 
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2.3. Density Matrices 

The behavior of the nuclear spin system is best described by an ensemble averaged 
density matrix, p{t) = \ilj(t)){ip(t)\, where li'it)) is the wavefunction of a single nuclear 
spin in the ensemble. The quantities measured in an NMR experiment, Mx{t) or Mz{t), 
describing the transverse or longitudinal nuclear magnetization, can be written in terms 
of the density matrix via the expression M^it) = (Ma) = Tr{Ma ■ pit)}- Details of the 
spectrum and relaxation rates are therefore contained in p{t). For a typical spin 1/2 
nucleus, the density matrix is given by: 

m = 

where T^^ and are the spin-lattice relaxation rate and the spin decoherence rates, 
respectively, p^'^ = exp{—li./kT)/Z is the density matrix in thermal equilibrium, and Z 
is the partition function. Clearly, the relaxation rates tend to bring the density matrix 
to its equilibrium, time-independent value. Spin-lattice relaxation is responsible for 
the time decay of the diagonal terms of the density matrix, whereas spin-decoherence 
relaxation is responsible for the decay of the off-diagonal terms. The spin-lattice 
relaxation rate provides the time-dependence of the component of the nuclear spin 
magnetization, whereas the decoherence rate T2~^ affects the transverse magnetization, 
M,{t). 

2.4- Relaxation Processes 

The spin lattice relaxation rate arises from fluctuations which couple nuclear spin levels 
which have a finite power spectral density at the nuclear Larmor frequency, ul = 'jHq. 
The most straightforward method to calculate T^^ is via first order time-dependent 
perturbation theory, which gives the Fermi golden rule: 

P^^f = ^\m^\f)Mco-Ef + E,) (9) 

where Pi^f is the rate of transitions between states \i) and |/), and Tii is the 
(time dependent) Hamiltonian giving rise to the spin lattice relaxation. In heavy 
fermion materials, typically Tii is given by hyperfine interaction (Eq. [7]), where the 
time dependence arises from the dynamical fluctuations of the spins S(rj,t), or from 
coupling to the orbital magnetization in a vortex lattice, which may fluctuate due to 
vortex motion. The principle of detailed balance assures that an ensemble of nuclear 
spins will reach equilibrium through a series of transitions between levels as described 
above. In equilibrium, the population of the nuclear spin levels can be described by 
a spin temperature, Tg, such that the population of each level. Pi, is described by 
the Boltzmann factor: Pj = exp{—Ei/ki,Ts)/Z, where Ei is the energy of the i^^ level, 
and Z is the partition function. One can show that the time dependence of the spin 
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temperature Tg = l/ksPs can be described as: 

^ = ^(/3.W-/5i.Q) (10) 

where Peq is the equihbrium spin temperature, and Tf ^ is the spin lattice relaxation 
rate [T6] . 

In metallic systems, the dominant coupling to the nuclei is the hyperfine interaction, 
and the spin lattice relaxation rate can be written in terms of the autocorrelation 
function of the dynamical hyperfine field, -ffhf (^) as: 

Tf 1 = 1- / {H,,,{t)Hy,,{0))e'^'^'dt, (11) 
^ Jo 

where c^o is the Larmor frequency. If we assume that H\i{{t) fluctuates randomly between 
±ifhf with autocorrelation function: {HM{t)Hhi{0)) = if^fC"*^'^, then we can write: 

This is a reasonable description of many processes involving phase transitions, where the 
fluctuation time scale changes by several orders of magnitude through a phase transition. 
The limit ujqt ^ 1 corresponds to motionally narrowed limit, as described below. Note 
that in this case, T{~^ is independent of frequency. 

2.4-1 ■ Block Equations By employing Eq. [8], one can derive expressions for the 
expectation value of the nuclear magnetization as a function of time. The expressions 
are the famous Bloch equations, given by: 

-^=7(MxH),+ (13) 
= 7(M X H), - — (14) 

^ = 7(M X H), - ^. (15) 

These equations not only make a connection to the classical behavior of a precessing 
magnetic moment in a magnetic field, but they reveal a more physical picture of the 
relaxation rates T^^ and Tg"^. Clearly, T^^ corresponds to a longitudinal relaxation 
of the magnetization along the z-axis, whereas T2 is a transverse relaxation in the xy- 
plane. Quantum mechanically, T2 corresponds to a decoherence time, in which the 
ensemble of precessing nuclear spins looses its coherence. Immediately after a pulse, 
one creates off-diagonal coherences in the density matrix. Over a time scale, T2, the 
nuclear wavefunctions are no longer in phase, and the precessing nuclear spins are no 
longer in step [I7j. 

2.5. Motional Narrowing 

When a nucleus experiences a static field, perhaps from an ordered magnetic phase, 
the NMR spectrum will reflect this internal field, usually as a broadening or a shift of 




Figure 1. Calculated spectra for a nucleus experiencing a dynamic hyperfine field ±/ihf 
with correlation time r. The graphs are offset vertically, with values of jhi^fT2 — 4, 
and T2/T — 100, 10, 8, 1, 0.1, 0.01 from the upper graph to the bottom graph. 



the resonance frequency in the absence of this field. On the other hand, if this field is 
fluctuating on a fast time scale, the effect will be averaged out and will not contribute 
to the NMR spectrum. This phenomenon is known as motional narrowing, and plays 
a key role in the NMR physics of many heavy electron systems. The relevant scale for 
the crossover from the static to the motionally narrowed limit is given by "jkhfT ~ 1, 
where hhf is the magnitude of the fluctuating field, and r is the correlation time of these 
fluctuations. This phenomenon is shown in Fig. [T], where the spectra are shown for a 
series of values of r, where the field hops randomly between ±/i/i/ with correlation time 
r. A detailed description of the general phenomenon can be found in [T6] . 

2. 6. NMR and magnetic phase transitions 

In a typical magnetic phase transition, both and the spectrum change dramatically. 
To illustrate this, let us assume that there is an internal magnetic field /ihf(^) from 
the magnetic order that is created by the hyperfine interaction at the nuclear site. 
Above the magnetic ordering temperature, hh{{t) fluctuates randomly from ±ho with 
an autocorrelation function (/^hf (^)^hf(O)) = /iqC"*/^. If the transition is second order, 
the temperature dependence of r is typically described by r(T) = ri(l — T/T^)"''^ for 
T < Tm and r(T) = T2{T/Tm — for T > T^. When critical opalescence sets in 

above the phase transition, T^"^ typically increases. This can be seen in Fig. [2], where 
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2 4 6 8 10 

Figure 2. The spin lattice relaxation rate versus {llJot)~^ showing the Bloembergen- 
Purcell-Pound (BPP) peak (Eq. [H]). 

T]~^ is plotted as a function of 1/r; i.e. high temperatures would correspond to large 
values of 1/r. As the correlation time grows and reaches a value on the order of the 
Larmor frequency, the fluctuating field becomes most efficient at inducing transitions 
between the nuclear levels and T^^ reaches a maximum. Fig. [3] shows the spectrum 
and the spin lattice relaxation rate as a function of temperature assuming t] = 1.1. 
Clearly, there are dramatic changes as the temperature goes through Tm, and static 
order sets in. This scenario is a good model for a second order antiferromagnetic phase 
transition, since the hyperfine field from the ordered spins is either positive or negative 
at the nuclear sites in the unit cell. However, in real systems the ordered moment grows 
continuously from zero, whereas in the case studied here it is finite and temperature 
independent. 

3. NMR in Metals 

3.1. RF Radiation in Metals 

NMR in metallic systems is experimentally challenging because the radiofrequency skin 
depth limits the effective volume of sample that can be measured. This skin depth 
affects both the excitation pulses (Hi) and the signal from the precessing nuclei within 
the bulk. The rf skin depth 6 = y^2pjujjj,, where p is the resistivity of the metal, is on 
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Figure 3. Spectra and relaxation rate for an ensemble of nuclei experiencing a dynamic 
hyperfine field ±ft,ht as a function of temperature. Here we consider a correlation time 
T = To{T/Tm — 1)^'' for T > Tm and rj = 1.1, where Tm represents a magnetic 
ordering temperature at which point the electron spins (and hence the hyperfine field 
at the nucleus) becomes static. Note that the temperature dependence of the spectra 
below Tm differ from typical materials at a second order transition, where the ordered 
moment (and the hyperfine field) grow continuously from zero. In this case, there is 
only a single value of | /ihf | . 

the order of a few microns for good metals at typical NMR frequencies, whereas sample 
sizes are typically on the order of a few mm. Nevertheless, the number of nuclei within 
the skin depth layer of a single crystal may reach the order of 10^'', which is sufficient 
to detect. Reducing the temperature and averaging several echoes (typically between 
10 to 10^) increase the signal-to-noise ratio (SNR) and make it possible to accurately 
measure spectra and relaxation rates. Powdered polycrystalline samples increase the 
surface to volume ratio and boost the SNR significantly, but can be complicated by 
orientation-dependent terms such as the quadrupolar interaction for nuclei with spin 
greater than / = |- To get around this, one can also work with aligned samples, in 
which the powder is mixed with epoxy and cured in a magnetic field. If the sample has 
an anisotropic magnetic susceptibility, the crystallites will align in the external field. 

3.2. Hyperfine coupling and Knight shifts 

In a Fermi liquid the dominant coupling is the contact hyperfine coupling to the 
quasiparticle spins, given by: 

^hyp = v4i ■ S 

= AhS, + ^Aii+S^ + LS+). (16) 

This coupling allows the nuclei probe both the static susceptibility, xo, and the 
dynamical susceptibility, x(q, a;), of the quasiparticles. The diagonal term AI^S^ gives 
rise to a static shift of the resonance frequency. The Knight shift, K, which measures 
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the percent shift of the resonance frequency from that of an isolated nucleus {uq = 'jHq) 
is given by: K = Axo/h'yfiB- In a Fermi liquid, Xo is given by the Pauli susceptibility, 
so K ^ AN{0) is temperature independent. 

This scenario works well for simple metals such as Li and Na, yet there are many 
cases where the Knight shift is more complex. In Pt, for example, there are multiple 
hyperfine couplings to the d- and sp- bands, and hence several contributions to the total 
shift [1^. In many-electron atoms, there is also a core-polarization term, in which the 
core s electrons acquire a population difference between the up- and down-spin states. 
This difference arises because the orthogonal eigenstates of the many-electron atom get 
mixed by the perturbing influence of the external field. In practice, it is difficult to 
estimate the contribution of a core-polarization term versus a purely contact term |19j . 
As a result, hyperfine couplings are usually taken to be material dependent parameters. 
In metals with local moments, such as rare-earth and d-electron systems, there is also 
a second hyperfine coupling to these moments that give rise to a strong temperature 
dependence of the total shift. 

3.3. Spin-flip scattering and spin-lattice relaxation 

The off diagonal terms of the hyperfine coupling can be written as: ^(/_|_S'+ -|- I^S^). 
This perturbation corresponds to a spin-flip exchange between the quasiparticle spin 
and the nuclear spin. These processes do not shift the resonance frequency, but do 
affect the dynamics of the nuclei. In fact, this process is crucial to bring about an 
equilibrium population distribution among the nuclear spin levels. In a Fermi liquid, 
as the quasiparticles scatter from one nucleus to another, they maintain essentially the 
same energy since the nuclear Zeeman energy is orders of magnitude lower than the 
Fermi level, and the quasiparticle Zeeman energy can be absorbed by states within fc^T 
of Ep- By using Fermi's Golden Rule one can show that the spin lattice relaxation rate 
can be written as: 



where the brackets indicate a thermal averaged correlation function. In a Fermi liquid, 
the states available for scattering the quasiparticles are those at the Fermi surface, and 
by a simple counting argument one can show that (JTTll can be written as: 



where Ef - Ei = hwo- Since /(S)(l - f{E)) ^ ksTdiE - Ep), we find that 
Tf ^ ~ Ty4^A^^(0). In other words, measurements of T^^ yield information about the 
square of the density of quasiparticle states at the Fermi level. Herein lies the power of 
NMR to probe the phenomenon of superconductivity: changes to the Fermi surface, such 
as the development of a gap, will be reflected in T^^. It is also immediately obvious that 
TiTK^ should be a constant in a Fermi liquid. This Korringa constant, JC = Tt'^h'-y'^ / ix\, 
is valid for non-interacting systems, and there are very few experimental systems for 




(17) 




(18) 
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which TiTK'^/K, is unity. The value of this ratio, in fact, is a measure of the strength 
of the quasiparticle interactions [20] . 

In 1957, Torn Moriya recognized that the expression ( |T71) can be rewritten also 
in terms of the dynamical susceptibility of the quasiparticles. By employing the 
fluctuation-dissipation theorem, Eq. (ITTj) can be written as: 



where the form factor A^(q) is the square of the Fourier transform of the hyperfine 
interaction, and the sum over q is over the first Brillouin zone. For a contact interaction, 
A^(q) is constant, but for more complex situations involving transferred couplings 
between neighboring sites, v4^(q) can have structure and may vanish at particular 
wavevectors. A q— dependent form factor can have profound consequences for the 
behavior of T^^ in materials. A notorious example is the difference in T^^ observed 
for the planar Cu and planar O in the cuprates [21]. Each nucleus has a different form 
factor and the dynamical susceptibility of this material is dominated by fluctuations at 
a particular wavevector, Q. Since v4^(q) vanishes for the O site but not for the Cu site, 
the two spin lattice relaxation rates have very different temperature dependences, even 
though they are coupled to the same degree of freedom, the Cu 3d^ S" = 1/2 spins. 

4. Probing the superconducting order via the hyperfine field 

4.1. Knight shift measurements in type-I and type-II superconductors 

The hyperfine interaction couples nuclear spins to quasiparticle spins, but in a 
superconductor at zero temperature, all of the 5* = | electrons are condensed as Cooper 
pairs in the superconducting condensate. The response of the nuclei depends critically 
on the spin of the Cooper pairs, which can be either even parity (spin singlet) pairing 
with net spin S* = or odd parity (spin triplet) with S = 1. For spin singlet pairing there 
is no coupling of nuclei to the Cooper pairs, and hence one expects all hyperfine fields to 
vanish at T = in the superconducting state. Therefore, measurements of the Knight 
shift, or the spin susceptibility, in the superconducting state offer powerful insight into 
the symmetry of the Cooper pairs. The spin susceptibility of a superconductor with 
spin-singlet Cooper pairs is given by the Yosida function: 



where Ns{E) is the density of states in the superconducting state, f{E) is the Fermi- 
Dirac distribution, and Xn is the (Pauli) susceptibility in the normal state. This function 
clearly vanishes as T — > 0, a result that is independent of the applied field direction. 

For spin-triplet pairing, the 5 = 1 Cooper pairs do couple to the nuclei, but the 
hyperfine field depends strongly on the orientation of the applied magnetic field with 
respect to the spin of the Cooper pairs. In other words, the behavior of the Knight 
shift below Tc depends on the projection of the S = 1 Cooper pair spin along H. To 




(19) 



q 
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NMR in heavy fermion superconductors 



14 



illustrate this anisotropy, we note that the superconducting order parameter can be 
written as: \^{k)) = A^{k)\ H) + Ao(k)(| Ti) + | Ti)) + A^(k)| TT)- Note here that 
\ipOi.)) spans the states {\Sz = +1), \Sz = 0), = —1)} in spin space. If we transform 
to a new basis {IS", = 0), \Sy = 0), \S, = 0)}, where ic = = 0) = (-| TT) + I ii))/V2, 
y = = 0) = (I TT) + I ii))/V2, and z = = 0) = (| Ti) + I iT))/V2, then we 
can write: |'?/'(k)) = dx{k)^ + dy{'k)y + dz(k)z. The order parameter is thus defined by 
a vector, d(k) = ''^^ ^^^'' x + ^^'^^.^^•' y + Aqz in spin space. For the case d x d* = 
(unitary state), d defines a direction in spin space given by the normal to a plane in 
which the spin expectation value remains finite. Typically, the orientation of the d- 
vector is pinned to the crystalline axes by the spin-orbit coupling. Depending on the 
components of d relative to the applied magnetic field, H, the expectation value of the 
spin either vanishes or remain finite. For example, in Sr2Ru04, the order parameter is 
commonly accepted to be d(k) = Aoz(A;^ ± iky), which corresponds to = +1, with 
the paired spins lying the the xy plane. In this case, if H lies along the z direction, 
will vanish since (Sz) = 0, but K^y will remain finite. Indeed, the Knight shift in the 
superconducting state of Sr2Ru04 shows exactly such behavior [22|. For further details 
on the triplet state in Sr2Ru04 the reader is referred to p3], and for a more general 
treatment of the order parameter in a triplet superconductor, see |24j . 

Measurement of the Knight shift in a superconductor is a non-trivial enterprise, 
and may seem contradictory at first glance. After all, how can one measure a resonance 
frequency in a material experiencing a Meissner effect? Furthermore, the magnetic 
susceptibility of a superconductor is dominated by the diamagnetic response of the 
orbital current response, and the spin susceptibility is only a small paramagnetic 
contribution. In fact, it is precisely because of the hyperfine coupling that one can 
extract the spin contribution in an NMR experiment, which cannot be done via bulk 
magnetization, or any other measurements. However, measurements of the Knight shift 
in a type-I superconductor are possible only in special situations in which the magnetic 
field can penetrate the sample, such as in the case of surface superconductivity. Some of 
the first NMR experiments on superconductors focused on Knight shift measurements 
of superconducting Hg colloids, where the surface area is very large [25]. In a type II 
superconductor, the situation is much clearer. In this case, the magnetization becomes 
position dependent in a vortex lattice, where M(r) = Morb(r) + Xspin-f^- Since the nuclei 
couple to the local internal magnetic induction, the resonance frequency will become 
position dependent as well: /(r) = ^{H + 47rM(r) -|- KH). In a type-I superconductor 
where M = —H/Att, the resonance frequency will be essentially zero except in special 
cases. However, in a type II superconductor, \M\ < H/An and so the resonance 
frequency is non-zero, although reduced from its value in the normal state. 

Experimentally, one often finds the resonance frequency shifts down in the 
superconducting state, as expected. The challenge is to extract the Knight shift from 
the somewhat complex spectra, since the dominant effect may be from the diamagnetic 
response of the orbital currents, rather than the suppression of spin susceptibility. 
The center of gravity of the distribution (Morb) = / ^orb(i")'^i' gives rise to a second 
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contribution to the shift of the resonance hne: -ft'orb = (^orb) / H, which is a complex 
function of field, the penetration depth, and the upper critical field ifc2- It is generally 
impossible to calculate the absolute value of (Morb)- Furthermore, there is also a 
demagnetization factor N/Air that shifts the net internal field inside the sample as a 
result of the diamagnetic magnetization, and suppresses the resonance frequency. This 
demagnetization field is a strong function of sample geometry, and depends on the 
direction of the external field with respect to the sample. The combined effect of K^rh 
and the demagnetization field can cause the resonance frequency to change dramatically 
in the superconducting state, even in the absence of a Knight shift (spin susceptibility) . 
In practice, then, it can be particularly challenging to determine what part of the 
measured shift arises from a suppression of spin susceptibility and what arises because 
of the diamagnetic response. However, in superconductors which have more than one 
NMR-active nucleus, one can sometimes decouple these two effects if one knows the 
hyperfine coupling to each [26] . A critical evaluation of any Knight shift measurement 
in the superconducting state needs to consider all of these effects in detail. A further 
challenge to measurements of the Knight shift at T = is the spin-orbit coupling. 
Spin-orbit coupling mixes spin-up and spin-down states such that the Cooper pairs are 
composed of a mixture of the two, and the ground-state spin susceptibility will approach 
a finite value at T = rather than vanish. The review by MacLaughlin discusses these 
effects in detail, as well as several experimental results [27] . 

Once the zero-temperature Knight shift is determined, and the orbital, 
demagnetization field, and spin-orbit effects have been evaluated, the remaining 
temperature dependent shift arises from the suppression of the spin susceptibility in 
the superconducting state. In principle, the detailed temperature dependence can be 
calculated with Eq. fl20|) . and compared with the density of states to extract quantities 
such as the magnitude of the superconducting gap, | A|, and the symmetry of the gap in 
k-space. In practice, however, the differences between the calculated shifts for different 
gap functions easily can fall within the noise. 

It should be noted that the SNR in the superconducting state is usually significantly 
reduced as a result of the enhanced screening of the rf excitation pulses. Whereas in the 
normal state, the rf skin depth is on the order of a few microns {B{r) ~ e-(^+^)^/''), in the 
superconducting state, the rf penetration depth is often only an order of magnitude less 
(S(r) ~ e-^'/^c-) [28]. For NMR Knight shift measurements, one requires single crystals 
rather than powders, and therefore this change of NMR-active volume can significantly 
reduce the signal, particularly just below T^. As the temperature is lowered below Tc, 
the Boltzmann factor will tend to enhance the signal. In a type-II superconductor, the 
rf pulses penetrate to within the Campbell depth, > A. In effect, the vortex lattice 
itself responds to the external time-dependent field, resulting in a greater rf penetration. 
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4-2. Spin-lattice relaxation in superconductors 

In a fully gapped spin-singlet superconductor at T = there are no fluctuations of 
Hhypit) to couple to the nuclei. At finite temperature, quasiparticle excitations couple 
to the nuclear spins, but since the excitations of a BCS state are actually superpositions 
of spin-up and spin-down quasiparticles, the hyperfine field experienced by the nuclei 
is different than in the normal state. This fact has profound consequences for the spin 
lattice relaxation rate in superconductors, and gives rise to the Hebel-Slichter coherence 
peak. To see this, we can write Eq. (fT7|) as: 

Tf ^ °^ y cl^c^^dkdk', (21) 

where c|^^ is the creation operator for the state |ks). This expression clearly describes 
a spin-flip scattering event weighted by the square of the hyperfine coupling. The 
Buglioubov excitations in a superconductor are coherent superpositions of the states 
|k 1) and | — k |), such that: 

"kl,2 = ^kCki ± l^-kC-ki, (22) 

where aj^^ 2 is the creation operator of a Bugoliubov excitation above the 
superconducting condensate, and Wk and are their respective amplitudes. Therefore, 
Tf^ can be rewritten as: 



T^^ (X y (wk'Mk + fk''yk)ttk'ittk2C?k(ik'. 



(23) 



This expression shows that the spin-lattice relaxation in the superconducting state is 
driven by spin-flip scattering of Bugoliubov excitations, but weighted by a coherence 
factor C+(k, k') = u^'Uk + fk'^k, a quantity that is determined by the superconducting 
condensate parameters. If the superconducting state were actually a simple two-fluid 
system rather than composed of coherent Cooper pairs, then the relative phase of fk 
with respect to would vanish, C+(k, k') would average to unity, and Eq. fl23l) would 
simply describe relaxation by a new fluid of al.^ quasiparticles. However, since the 
Bugoliubov excitations are indeed coherent superpositions of the simultaneous spin-up 
and spin-down electrons, the coherence factor is greater than unity, which has profound 
consequences for the relaxation rate, as shown below. 

Eq. (I23l) can be simplified by noting that the number of Bugoliubov quasiparticles is 
given by the Fermi-Dirac distribution, and employing the superconducting state density 
of states: 

Tf^ oc C+{E,E')NBcs{E)fiE)NBcsiE')il - f{E'))dE, (24) 

where the density of states of an s-wave (L = symmetry of the Cooper pairs) is given 
by: 

iVBcs(i?) = -^=^==N^{E) (25) 
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and the coherence factor is: 



EE' 



(26) 



2A2' 



Since N^csiE) has a singularity at = A and C+{E,E') is non-zero, expression 
leads to a singularity of T^^at T = Tc. This is the Hebel-Slichter coherence peak, first 
observed experimentally in superconducting Al as an unexpected peak in Tf^ just below 
Tc [2]. The original data measured by Masuda and Redfield is shown in Fig. HI which 
clearly shows the dramatic increase in T^^ just below Tc. The intensity of the peak 
is limited in practice by both a smearing of the density of states as well as spin-orbit 
scattering effects. The spin-orbit interaction mixes spin-up and spin-down states with 
different |k) states, so the Bugoliubov states in Eq. [22] are modified to include a sum 
over k states, as well as the coherence factor C+(k, k'). In practice, the coherence factor 
can be rather difficult to observe experimentally, particularly in materials with heavier 
elements where the spin-orbit interaction is larger. Fig. ([5]) shows data for a number 
of s-wave superconductors, and clearly shows how the intensity of the coherence peak 
can be strongly suppressed, particularly in larger fields. This suppression may be due 
in part to localized excitations within the vortex cores [29]. The reader is referred to 
[27] for further details. It is important to note that this peak is relatively hard to see in 
most materials and that its absence is not necessarily an indication of unconventional 
superconductivity. 

At this point, it is worth noting that the acquisition of these data during the 
time period 1957 - 1961 represents an experimental tour de force, since Al is a type-I 
superconductor with a transition temperature of 1.17 K. Although other nuclei were 
available with higher T^'s, or even ^^^In with an NQR frequency of 7.54 MHz and 
Tc = 3.41 K, Al was chosen because the low Z reduced potential effects of spin-orbit 
coupling (In NQR was performed ten years later by Butterworth and MacLaughlin 
[30]). Hebel and Slichter at the University of Illinois and Masuda and Redfield at 
IBM both tackled the problem independently by constructing an ingenious field cycling 
technique to let the nuclei relax in zero external field (in the superconducting state), 
while measuring the nuclear magnetization in an applied field in the normal state [SIE]. 
Hebel and Slichter were able to reach temperatures ~ 0.94 K by pumping on liquid He-4. 
Masuda and Redfield reached 0.35 K by pumping on He-3 using a special He-3 cryostat 
that had just been developed [31]. Aside from the rigorous and technically challenging 
field and temperature requirements, even simple NMR measurements were challenging 
by today's standards, where standard NMR spectrometers consist of computer controlled 
digital radio transmitters and receivers. Indeed, many frequency measurements were 
made by using a Lecher-wire assembly, in which radiofrequencies were measured by 
setting up standing waves on parallel wires, and the nodes of the standing waves were 
measured by connecting a lightbulb across the pair at various distances. 

For T < Tc/3, Eq. ( 12^ reduces to ~ exp(— A/Zc^T). In other words, the spin 
lattice relaxation rate exhibits an activated temperature dependence, as expected in the 
presence of a uniform gap function. This behavior, however, is only present if the gap is 
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Figure 4. The relaxation rate of tlie ^^Al in the normal and superconducting states, 
reproduced from [3]. Note the dramatic coherence peak in T^^ just below Tc- The solid 
red line is a fit to T^^ ^ exp(— A/ZcsT) as described in the text, with 2A = SAksTc, 
and the dashed green line is a plot of T^^ ^ T^, as expected for a nodal superconductor. 



isotropic in k-space, i.e., for s-wave superconductivity. If L > and Ak develops nodes, 
then the density of states A^bcs(-E') is modified. For hne nodes in an L = 2 (d-wave) 
superconductor, Nbcs{E) ~ E/Aq for E <^ Aq, where Aq is the maximum gap value. 
In this case. Eg. (1241) leads to T^^ ~ T'^ for T < T^/S. For point nodes, as may be 
the case for L = 1 (triplet) pairing, then T^^ ~ |24]. These dramatic differences 
in temperature dependence arises because there are always quasiparticle excitations 
available at the gap nodes. Measurements of the temperature dependence of T^^ can 
help to distinguish between these different types of pairing. It should be noted, however, 
that Tf^ is only sensitive to the magnitude of of the gap, not the sign. In principle, 
if the gap is anisotropic s-wave, then the density of states can mimic that of a true 
superconductor with nodes in the gap. In the case of the cuprates, the quantity T2G, the 
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Figure 5. The relaxation rates of tiie skutterudite compound LaFe4Pi2 [32j, for 
MgB2 [33 , elemental Al [3], and elemental In [30] . These materials are all s-wave 
superconductors, yet the visibility of the coherence peak varies dramatically. 



Gaussian component of the spin-spin relaxation rate, is sensitive to the change of sign of 
the gap function in k-space. T2G is a measure of the indirect couphng between nuclei as 
mediated by the real part of the electronic spin susceptibility, x'i^^ 00 = 0). This quantity 
is modified below Tc and is a function of Ak±q. Since T2G involves the coupling of two 
different nuclei to x'(q, c<j = 0), it can distinguish between an anisotropic s-wave and a 
true d-wave superconductor [M]. Indeed, measurements in YBa2Cu408 are consistent 
with calculations based on d-wave, rather than anisotropic s-wave, superconductivity 
[35] . This connection between T2G and x'(q, u; = 0) is a special case that happens to 
apply to the Cu nuclei in the cuprates. The hyperfine couplings of the Cu are such 
that the indirect nuclear-nuclear coupling is Ising-like with a particular direction (the 
c direction for the Cu). In this case, the form of the echo decay is Gaussian, and is 
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characterized by the quantity T2G, which can be directly related to x'(q, c<j = 0) [55] . 
In general, though, the indirect couplings are weaker and more isotropic than in the 
cuprates, and consequently there is little information to be extracted from the form of 
the echo decay. 



5. Probing the vortex lattice 
5.1. Field distributions 

In a type II superconductor, the inhomogeneous distribution of supercurrent in a vortex 
lattice gives rise to a periodic array of local orbital magnetization, Morb(r). The orbital 
magnetization in a vortex lattice couples to the nuclear spins via a chemical shift: 
'Horh = 47rl ■ Morb(r). Typically a chemical shift is given by an on-site diamagnetic 
orbital susceptibility in a solid or molecule. In a superconductor, the diamagnetic 
response of the supercurrents is a macroscopic phenomenon, with an inhomogeneity 
set by the scale of the magnetic field. The unit cell of a triangular vortex lattice is given 
by: rvL ^ ^/^o/Hq ~ ^jHc2/ H^, where ^ is the coherence length. Since Tvl 3> a, 
the unit cell length, the local orbital shift (and hence the NMR frequency) varies on a 
length scale much larger than the vortex lattice. Each nucleus resonates at the local 
magnetic field /(r) = 7|ifo + 4'/rMorb(r)|, so the NMR spectrum consists of a histogram 
of local fields in the vortex lattice, V{f) ~ l/|V/(r)| (see Fig. [6]). This characteristic 
spectrum, the "Redfield pattern", has been used extensively to determine properties of 
the vortex lattice in a number of superconductors [371 IMl l2Z] • Ideally, one can measure 
the spectrum V{f) and fit the data to an analytic expression, but this is not realistic 
in practice. Typically other inhomogeneities limit the spectral resolution such as the 
magnet inhomogeneity. Fig. [7] shows the ^^P spectrum in the mixed state of LaRu4Pi2 
[39] . Although the spectrum is broadened, a clear tail is evident at the upper frequency, 
corresponding the the classic Redfield pattern. 

As discussed above, there can be errors in determining the absolute value of the 
shift, given by the first moment of the spectrum, because there are several quantities 
that affect the resonance frequency in the superconducting state. The linewidth, or the 
rms second moment of the spectrum, is often more useful and less prone to extrinsic 
effects. The second moment is given by 

((//7 - Hof) ^ 0.00371^1 (27) 



for H <^ Hc2 and 



for H close to ifc2, where A is the penetration depth, and k = A/^ is the 
Ginsburg-Landau ratio. Therefore, measurements of the lineshape can provide 
direct information about the superconducting penetration depth (and hence the 
superconducting condensate density Ug ~ 1/A^) as a function of temperature [26], [38]. 
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Figure 6. The field distribution 'P{h) versus local field h = H — Hq in a series of 
external fields Hq — 0.5, 1.0, 5.0, 10.0 and 15.0 T for a superconductor witli Hc2 = 36.0 
T and penetration depth A — 2700 A. The different fields have been offset vertically 
(with increasing field) for clarity. 



Vortex matter can undergo phase transitions, and NMR can probe the physics of 
the vortex lattice phase diagram. In two-dimensional layered systems such as the high 
temperature superconductors, the pancake vortices can be either coupled or decoupled 
between layers, and undergo a vortex liquid to solid transition [IQ]. Since the vortices 
can move in the liquid state, they create a dynamic field 5M(r, t) at the nuclear site 
that can lead to relaxation. If the autocorrelation time of the longitudinal fiuctuations, 
r, defined as {SMz{r,t)5Mz{r,0)) ~ e~*/^, is short compared to l/'y5Mz, then the 
spectrum will be motionally narrowed, and will not show the characteristic Redfield 
pattern. When the vortex lattice cools and undergos a liquid to solid transition, the 
correlation time r will become much longer, in which case the condition I'^SM^t ^ 1 
will be satisfied. In this case, the spectrum will no longer be motionally narrowed, and 
the Redfield pattern will emerge. This technique was applied to map the vortex lattice 
phase diagram in YBa2Cu307 [H]. Also, in many instances anomalies have been found 
in Tj~^ and Tg"^ at characteristic temperatures inside the superconducting state, which 
are probably related to motions of the vortex lattice [i2l [35] . A solid vortex lattice can 
also experience overdamped vibrational modes as a results of the finite moduli in the 
lattice [m HO] . These modes can give rise to fiuctuating fields at the nuclear sites close 
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to the vortex cores, and may contribute a significant portion to the local spin lattice 
relaxation rate. For this reason, it is advantageous to measure Tf ^ at zero field, in 
the absence of vortices. Typically this is done in an NQR experiment, if a quadrupolar 
nucleus is available and has a reasonably large quadrupolar frequency (> 1 MHz). 

5.2. Interaction with radiofrequency irradiation 

The radiofrequency NMR pulses used for exciting the resonant nuclei can also drive 
motion of the vortex lattice. In order to investigate the effects of the rf pulses on 
the vortex lattice itself, researchers at UCLA have made systematic studies of how 
the rf pulses can anneal a distorted vortex lattice [Mj. In these studies, a type-II 
superconducting material is field-cooled, then rotated by an angle on the order of a few 
degrees with respect to the static field. This rotation strains the vortex lattice, giving 
rise to a large inhomogeneous field distribution. In a standard spin-echo experiment, 
the field inhomogeneity is refocused, and the original signal is recovered at a time 2r 
from the first pulse. In the case of the strained vortex lattice, the magnitude of the 
echo is reduced from that of the unstrained lattice. The reason is the refocusing tt 
pulse partially anneals the vortex lattice, and changes the local field at points in the 
vortex lattice where the positions of the vortex lines have changed. In these regions, 
the signal is not refocused. However, after a series of pulses, the full signal intensity 
can be recovered. This experiment suggests that the rf pulses used in NMR of type II 
superconductors do not significantly affect the equihbrium field distribution created by 
field-cooled conditions. 

The rf pulse annealing experiments also suggest that significant field distributions 
can exist when the vortex lattice is strained. Therefore, NMR experiments that are 
designed to probe a field distribution (such as measurements of the penetration depth 
described above) should always be conducted under field-cooled conditions. However, 
in order to extract an NMR spectrum, it is often easier and more straightforward to 
take data by sweeping frequency at constant field rather than take data by sweeping 
frequency at constant field. Obviously the latter is going to be more accurate. Indeed, 
the former technique is likely to reveal extrinsic details of the field inhomogeneity rather 
than intrinsic details of the physics of the vortex lattice. If the field is swept inside the 
Abrikosov vortex lattice state, then the vortices will either enter or depart the sample 
from the ends, and the magnetization will be inhomogeneous on a macroscopic scale 
with a magnitude determined by the critical current density, Jc, as described in the Bean 
model [15] . In this case, it is not straightforward to determine if quantities such as the 
second moment of the NMR line are probing intrinsic properties such as the penetration 
depth, or extrinsic properties that depend on the sample geometry. Unfortunately this 
phenomenon is not widely recognized or discussed in the NMR community, and non- 
specialists are cautioned that reports of these types of experiments should always include 
a careful discussion of the experimental conditions. 



NMR in heavy fermion superconductors 



23 



1.4 



1.2 - 



1.0 - 



8 0.8 



0.6 
0.4 



0.2 - 



0.0 




16.88 16.92 16.96 

Frequency (MHz) 



17.00 



Figure 7. The '^^P spectrum (solid line) and the spin lattice relaxation rate (•) in the 
mixed state of superconducting LaRu4Pi2 {Tc — 7.2 K) at 2 K and 5.2 T. Reproduced 
from m. 



5. 3. Imaging experiments in vortex lattice field gradients 

A significant advantage of NMR experiments in a vortex lattice is tlie possibility 
of real-space imaging because of the large field gradients intrinsically present in a 
vortex lattice. Since the resonance frequency in the Redfield pattern corresponds to 
a particular location in real space relative to vortex cores, one can selectively measure 
the relaxation rates at different parts of the vortex lattice [IHIHTIHH]. This technique is 
particularly enlightening in the case of d-wave superconductors, where the local density 
of quasiparticle states varies spatially because of the so-called "Doppler" shift. The 
energy of a quasiparticle excitation which is moving as part of a supercurrent will have 
an energy: 

^k= + A^ + ^v,-k, (29) 
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where is the normal state dispersion, and is the local supercurrent velocity. The 
second "Doppler shift" term was first recognized by Volovik [19], and is important 
because it gives rise to an enhancement of the local density of states. At T = the 
point nodes of a nodal superconductor become finite arcs in k— space because of this 
Doppler shift. Since the supercurrent velocity varies spatially in a vortex lattice with 
a periodicity commensurate with the local field, there is a correspondence between the 
local resonance frequency and the local density of quasiparticle states, and hence the spin 
lattice relaxation rate. Several groups have taken advantage of this result to measure 
the properties of the vortex lattice in the cuprates [IHl HTl HH [501 EB 129] • These 
experiments generally find a strong spatial dependence of Tf^ outside of the vortex 
cores, as predicted. At short distances, essentially within the vortex core, T{'^ exhibits 
very different behavior in the case of YBa2Cu307 in high fields, and in Tl2Ba2Cu06+5 
compounds. This behavior inside the vortex cores (within a radius ^ of the center of the 
vortices, where ^ is the coherence length) is suggestive of antiferromagnetic fluctuations, 
and has been taken as evidence for antiferromagnetic vortex cores in the cuprates. 
Subsequent inelastic neutron scattering data seem to support this scenario [52], [53]. The 
coexistence of two order parameters is possible since the superconducting order is locally 
suppressed within the vortex cores. A similar phenomenon is probably at play in the 
heavy fermion superconductor CeCoIus [51]. 

Measurements of the local density of states in an s-wave superconductor can also 
reveal strong local spatial variations. Recently, Ishida et al. measured the Redfield 
pattern in the skutterudite superconductor LaRu4Pi2 (Tc =7.2 K), which is an s-wave 
superconductor. Fig. [7] shows a strong frequency dependence of Tf ^ similar to that 
observed in d-wave superconductors [46]. Clearly, the strong spatial dependence in this 
case cannot be attributed to the Doppler shift fl29l) . but must arise either from localized 
bound states in the vortex cores, or from vortex motion. 

6. NMR in the Heavy Electron Materials 

Superconductivity in the heavy fermion materials tends to emerge only under certain 
structures, and in particular regions of the phase diagram. In the case of the heavy 
fermions, much of the relevant physics can be captured by the Kondo lattice model 
in which a lattice of f-electron moments is embedded into a background of conduction 
electrons (Fig. [8^) [55]. The nature of the ground state is highly sensitive to the 
magnitude of the Kondo exchange interaction, J, and the conduction electron density 
of states, A^(0). When J and A^(0) are small, the f-moments can polarize the conduction 
electron medium, leading to an effective RKKY interaction between the moments and 
hence a magnetically ordered ground state. On the other hand, when J and N{0) are 
large, the conduction electrons form singlets with the f-moments, reducing the effective 
moment and leading to a spin-liquid ground state with no long range order. Between 
these two extremes there is a quantum phase transition between an ordered and a 
disordered state (Fig. [Hb). New broken symmetries often emerge in the vicinity of this 
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Figure 8. (a) In the Kondo lattice model a lattice of nearly localized spins (5/) 
interact with a sea of conduction electrons (Sc), via a Kondo coupling (J), (b) The 
phase diagram of the Kondo lattice, T versus JiV(O), where A^(0) is the density of 
conduction electron states. 



Quantum phase transitions have been focus of much theoretical and experimental 
interest in the past decade [56l [571 EH [Ml [60]. The critical fluctuations of a quantum 
phase transition are driven by quantum rather than thermal effects, and the observed 
behavior of many bulk observables near a QCP deviates strongly from conventional 
theories of phase transitions [611 152] • The Kondo lattice QCP described above is 
responsible for much of the new physics and emergent behavior exhibited by heavy 
fermion compounds. In these systems, the Kondo interaction, J oc V"^, is a strong 
function of the hybridization, V, between the conduction electron and /-electron 
wavef unctions. This hybridization can be modified by hydrostatic pressure or chemical 
doping and hence the ground state can be tuned over a broad range of phase space 
(Fig. [8]d). Consequently, heavy fermion materials offer an unprecedented testing ground 
to investigate the emergence of non-Fermi liquid behavior and the properties of novel 
ground states in the vicinity of a QCP. Although these systems are metallic, observed 
bulk properties such as the resistivity and specific heat disagree with the predictions of 
Fermi liquid theory, which is normally quite robust. This breakdown of Fermi liquid 
theory represents a significant challenge to theory and has been the subject of extensive 
experimental work in recent years [591 160] . 

NMR measurements can play a central role in probing the low energy degrees of 
freedom that emerge in the vicinity of a QCP. As the ground state evolves from an 
ordered magnetic state to a disordered state, the excitations above this ground state 
should change dramatically. These changes should be reflected in the low energy spin 
fluctuations, and hence in quantities such as Tj~^ Indeed, different theoretical models 
for quantum phase transitions can make specific predictions for the behavior of T^^ 
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Figure 9. Tj^^ versus temperature at various pressures in CeCoIus, reproduced from 
[66] . The strong pressure dependence may reflect the evolution of the low energy spin 
degrees of freedom as the system is tuned away from a QCP near P = GPa. 

[63] . The theory of local quantum criticality developed by Si and coworkers explains the 
unusual dynamical susceptibility, ^(q, ^, T) measured in some compounds, and predicts 
T]~^~ constant for sufficiently low temperatures in two dimensions |56]. The subsequent 
observation of exactly this behavior by Si NMR in YbRh2Si2 has helped identify this 
compound as a prime candidate for this theory [6^ . 

On the other hand, the situation is not so clear in CeCoIns, which is a particularly 
important material because it is a superconductor, and lies close to a QCP at ambient 
pressure [65]. Indeed, NMR measurements of the spin lattice relaxation rate as a 
function of pressure clearly show dramatic changes in the low energy degrees of freedom 
that are expected in the vicinity of a QCP (see Fig. [9]) [66]. The unusual T^/^ behavior 
observed at ambient pressure has been seen as evidence of quantum critical fluctuations 
[63| [66] . Although the strong pressure dependence must reflect the evolution away 
from a QCP, it is not clear whether the changes could be related to changes in critical 
fluctuations that persist up to 100 K, or perhaps changes in the hyperfine coupling itself. 
This material is just one of many examples in which the behavior in the normal state 
presents several challenges. 

6.1. The Knight shift anomaly 

In the normal state of heavy fermion materials, the magnetic susceptibility x usually 
increases strongly with decreasing temperature. In many instances, the Knight shift 
K fails to track the temperature dependence of x below a crossover temperature 
T* ~ 10 — 100 K (see Fig. [TOl) . This Knight shift anomaly is ubiquitous among the 
heavy fermions, yet is not fully understood, and suggests a priori that the hyperfine 
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coupling itself may be temperature dependent. This scenario would mean that all of 
the framework for understanding the Knight shift and Tf^ in the superconducting state 
(and indeed the normal state) would not be robust, leaving one with little confidence 
that the NMR data can be interpreted in any meaningful way. 

There are three theories of the Knight shift anomaly that have been proposed in 
the literature: two of these scenarios rely on a more local picture, and imply that the 
hyperfine coupling itself changes. The first, postulated by Cox and coworkers shortly 
after the discovery of the Knight shift anomaly in CeSns [67j, assumes that the f- 
electron spins are screened at temperatures T < T^, leading to a suppression of the 
bulk susceptibility, x [68j. The Knight shift, on the other hand, measures the local 
susceptibility and does not sense the screened moment. Although this scenario could 
explain the CeSns data, it fails to capture the behavior of the anomaly in several other 
heavy fermion systems where K x below T*. 

Another scenario was originally postulated by Yasuoka and Fisk to explain the 
Knight shift data in CeCu2Si2, in which the hyperfine coupling between the ligand 
Cu and Si nuclei and the Ce 4f electron spin depended on the crystal field states of 
the J = 5/2 multiplet [69j. As the temperature varies, the population of the CEF 
levels change and hence the relationship between K and x is modified. In this case, 
the local hyperfine field at a particular nucleus will depend the particular electronic 
configuration on the coupled 4f atom, which fluctuates on a fast time scale. The effect 
over the ensemble of nuclei would be a motional narrowing of the line, with an effective 
average field weighted by different hyperfine couplings to the susceptibilities of each of 
the CEF-split doublets. Data in CeCoIus was originally analyzed in a similar fashion 
[70] . but was ruled out later because T* = 50 K in this compound does not correspond 
to any of the CEF splittings, but rather coincides with the coherence onset temperature 
where many of the other observed quantities, such as the resistivity, change character. 

To date, the scenario that is most successful in capturing the behavior of a broad 
range of materials is the two-component hyperfine scenario, which asserts that the 
anomaly arises due to coherence in the Kondo lattice. In the Kondo lattice model, 
there are two sets of electron spins that contribute to the susceptibility: the conduction 
electron spins, Sc = 1/2, and the local moments of the f-electron, Jj. Typically, the 
f-electrons are localized at high temperature, giving rise to strong Curie- Weiss behavior. 
However, these conduction and f-electrons experience several interactions which modify 
their behavior at low temperature. The localized /-electron Jf multiplet is split by a 
crystalline electric field (on the order of T^cef ~10 meV), so that for T ^ T^cef/^b 
the magnetic behavior of this degree of freedom can be treated as a pseudospin doublet 
with an effective spin Sj = 1/2 and a g value that depends on the details of the 
crystal field parameters. Secondly, the local moments and conduction electrons are 
coupled via a Kondo interaction T^Kondo ~ J^c ■ Sf. This interaction gives rise to 
two different effects: Kondo screening of the local moments, and exchange interactions 
(Rudemann-Kittel-Kasuya-Yosida, RKKY) between the local moments. The Kondo 
interaction leads to a characteristic temperature Tx = N {fi) / k^e'^^ '^^^'^\ where A^(0) is 
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the density of states at the Fermi leveL Below this temperature, the f spins are screened 
by the conduction electron spins, reducing their effective moment. On the other hand, 
the Kondo interaction also gives rise to an indirect exchange between local moments 
via polarization of the conduction electron medium. The characteristic temperature 
for the RKKY interaction is given by Trkky ~ {•J^{^)Y- These two effects lead to 
two very different ground states: a spin liquid when the Kondo screening dominates, 
and long range magnetic order when the RKKY interaction dominates. Several years 
ago, Doniach considered the phase diagram of such a lattice of f-electrons, the Kondo 
lattice, and found that below a critical value of JN{0), the system would undergo long 
range magnetic order of the f-electrons, whereas for greater values the system remains 
disordered jTT]. This simple model describes some the essential physics, but fails to 
capture many details particularly in the vicinity of the quantum critical point (QCP), 
where superconductivity can emerge and the behavior of the normal state cannot be 
described as a Fermi liquid [57] . Despite many years of intensive theoretical work, there 
still remains no complete theory for the behavior of these heavy fermions in the vicinity 
of the QCP. 

Recently, Pines and coworkers have taken the approach that the essential physics 
of the Kondo lattice can be elucidated by characterizing the various materials in terms 
of a common phenomenology [721 [TSj [71|. This two-fluid picture assumes that below a 
temperature T*, the Kondo lattice is described by a coexistence of both local moment 
and heavy electron behavior. T* is a material dependent temperature on the order of 
10 - 100 K. This scenario offers an appealing microscopic picture of the Knight shift 
anomaly. For T > T*, K is linearly proportional to x- Since the dominant contribution 
to the susceptibility is from the local moments, Sf, at these temperatures, we conclude 
that the largest hyperfine field arises from the local moments. In general, however, there 
can be also an on-site hyperfine interaction to conduction electrons. A more complete 
description of the hyperfine interactions is given by: 



where A is an on-site hyperfine tensor interaction to the conduction electron spin, and B 
is a transferred hyperfine tensor to the f spins [15]. Note that we consider here nuclear 
spins on the ligand sites, i.e., not on the f atom nucleus. A similar scenario is present in 
the high temperature superconductors [13J. It is likely that the mechanism of the on-site 
coupling A is probably due to a combination of core polarization, and polarization of 
unfilled orbitals with s and p symmetry. The mechanism of the transferred interaction, 
B may arise as a combination of orbital overlap between the f electron states and the 
ligand atom wavefunctions, and either a core polarization of the ligand atom or a dipolar 
interaction between the ligand nucleus and the polarized ligand atom p- and d- electrons. 
For the remainder of this article, we assume that the hyperfine parameters are material 
dependent constants, and will not discuss their microscopic mechanism further. 

Given the two spin species, Sc and Sf, there are three different spin susceptibilities: 




(30) 
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Xcc = (ScSc), Xcf = (ScSf), and Xff = {^fSf). The full expression for the Knight shift 
is given by: 

K{T) = AxcciT) +Ia+Y.bA Xcf{T) + Baff{T). (31) 

\ i&iri / i&m 

where we have absorbed the g-factors into the definition of the hyperfine constants and 
dropped the tensor notation for notational simplicity [I5]. The bulk susceptibility is 
given by: 

X(T) = Xcc{T) + 2xc/(T) + Xff{T). (32) 

Note that ii A = Ylienn^i^ then oc % for all temperatures. In general, this relation 
does not hold, and if Xcc(T), Xcc(T), and Xcc(T) have different temperature dependences, 
then the Knight shift will not be proportional to susceptibility, leading to a Knight shift 
anomaly at a temperature T*. 

Since a complete solution of the Kondo lattice remains a challenge for theory at 
present, there is no description of the temperature dependence of Xcc(T), Xcf(T), and 
Xff{T). However, using empirical observations Nakatsuji Pines and Fisk, and later 
Yang and Pines, developed a phenomenological two-fluid picture of the susceptibility 
and specific heat in CeCoIns which also works well for understanding the Knight shift 
anomaly [72l [151 [751 [73] . In this picture, Xcf{T) ~ (1 - T/T*) log(T7T) for T < T*, 
and Xcf{T) ~ for T > T*, and we assume Xcc is negligible for all temperatures. This 
relation appears to hold for all heavy fermion and mixed valent compounds measured 
to date down to the relevant ordering temperatures; in fact, Xcf scales with T/T* 
in all instances. For a detailed discussion we refer the reader to [IHl [76]. Recent 
theoretical work and DMFT calculations lend support for the two- fluid picture [771 [78]. 
but a complete microscopic description is still lacking. Clearly, in the case of the 
heavy fermions it is the heavy quasiparticles that form the superconducting condensate, 
and the hyperfine coupling to these heavy quasiparticles is modified by the onset of 
coherence. Any experiment measuring the NMR properties in the superconducting state 
of a heavy fermion system must therefore somehow address the normal state Knight shift 
anomaly at some level. 

6.2. Superconducting order parameters 

6.2.1. CeCuQ,Si2 Since the discovery of the first heavy fermion superconductor, 
CeCu2Si2 in 1979, NMR has played a crucial role in determining the symmetry of 
the superconducting order parameter. The Knight shift of the Cu in this system is 
suppressed below T^., for fields oriented both parallel and perpendicular to the tetragonal 
c-axis (see Fig. [TT]) [79]. Note, however, that for the field along the c-direction, the 
shift increases below T^. This suggests that somehow the spin susceptibility increases 
below Tc, contrary to the behavior expected in a spin singlet superconductor. In fact, 
in the normal state between T^ and T* ~ 170 K, K^. decreases whereas Xc increases 
with decreasing temperature. Naively, one might interpret this as a negative hyperfine 
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Figure 10. The Knight shift of the In(l) (red squares) and the Co (blue circles) in 
CeCoIns, compared with the bulk susceptibility, x{T) (solid line), in the c— direction. 
INSET: Kc{ln{l)) and Kc{Go) versus Xc with temperature as an implicit parameter. 
The solid lines are linear fits to the data for T > T*. T*, the temperature where K 
and X diverge, is approximately 50 K in this case. 

coupling, so that below Tc when Xspin decreases (and vanishes in a spin singlet), the 
net effect on the total shift increases. A more sophisticated interpretation in the two 
component model is that {A — B)xcf{T) is negative in the normal state, and that \xcf\ 
decreases below Tc. However, it is not straightforward to associate Xcf with {ScSj) in 
this situation, as the character of the spins to which the nuclei are coupled in this heavy 
electron state are not simply described by Sc and Sj, but rather in terms of quasiparticle 
excitations of a highly correlated electronic state. Clearly, the susceptibility of these 
degrees of freedom is suppressed below Tc, but how the character of the local moments 
changes in the superconducting state is poorly understood. Spin lattice relaxation 
measurements in the superconducting state of CeCu2Si2 reveal a power law dependence; 
thus combined with evidence for spin-singlet pairing, it is likely that this material is a 
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Figure 11. The Cu Knight shift in CeCu2Si2. The data are reproduced from [79ll69j. 
The increase in Kf, below is consistent with the behavior above Tp, similar to that 
observed in CeCoIns. 



d-wave superconductor [80] . 



6.2.2. UPt^ - spin triplet superconductivity UPts is a highly unusual superconductor 
that exhibits multiple phases within the H — T phase diagram, that are most likely 
associated with different symmetries of the superconducting order parameter (see inset 
of Fig. [12]), and the majority of the data suggest that this is a spin-triplet superconductor 
[8T1[82]. The multiple phases possible below strongly suggest an order parameter with 
internal degrees of freedom. As described above, the response of the spin susceptibility 
below Tc depends on the direction of d relative to the field Hq. In other words, the spin 
susceptibility of the superconducting condensate must be described by a tensor that may 
have off-diagonal components. In a single crystal of UPta, Knight shift measurements 
of the Pt clearly showed no change between the normal and superconducting states for 
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Figure 12. The Knight shift of the Pt in UPta, reproduced from [83] ■ Open points 
correspond to the superconducting state, whereas soUd points are in the normal state. 
The Inset shows the H — T superconducting phase diagram for field H along the c 
direction. 



both field orientations [83]. This result is clearly incompatible with spin singlet pairing, 
and furthermore suggests that the spin polarization of the condensate is not pinned to 
the lattice, as was found in Sr2Ru04 [22], but rather follows the direction of the applied 
field, maintaining the same anisotropy as the normal state. 

6.2.3. CeCoIn^ and non-Fermi liquid behavior CeCoIns is particularly interesting 
heavy fermion superconductor because it shows a rich behavior of phenomena associated 
with the proximity of the system to an antiferromagnetic quantum critical point. The 
superconducting state phase diagram of this material exhibits an unusual field-induced 
transition to a new thermodynamic phase that lives only within a narrow regime between 
a field H*{T) and Hc2. This field induced phase will be discussed in detail below. For 
fields H < H*(T), K(T) and T{~^ behave in a similar manner to those in CeCu2Si2. The 
Knight shift of the In(l), In(2) and Co all change below Tc, but the sign of the change 
varies from site to site [70l|8l]. The In(l) shifts down below in both orientations, the 
Co and In(2) shift up for fields along c and down for fields along ab, and the ln{2)± shifts 
down for fields along ab (the In(2) splits into two distinct sites, In(2)|| and In(2)_L for 
fields along ab because it is in a low symmetry site). Empirically, if K{T) is increasing 
above Tc, then it decreases below, and if K{T) is decreasing above Tc, then it increases 
above Tc. As in the CeCu2Si2 case, this behavior can be understood in terms of a 
negative hyperfine coupling, and hence the spin susceptibihty decreases below Tc. 
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Figure 13. The Knight shift of the In(l) and the In(2) in CeCoIns, reproduced from 
[70 j . Note that the shift of the In(l) increases below Tc, similar to that of the Cu in 
CeCuzSia. 



The detailed behavior of the Knight shift as a function of field in the 
superconducting state of this material shows a dramatic increase with increasing field 
|85j . Xs{T = 0,H) should vary linearly with field for a d-wave superconductor because 
the Zeeman shift of the quasiparticle dispersion leads to a linear increase of the density 
of states at Ep [86]. In other words, the point nodes develop into arcs in k-space, as 
in the case of the Doppler shift. K{T = 0, H) of the In(l) shows a gradual increase, 
and eventually gets quite steep as H approaches H*{0) and Hc2 in this material. This 
behavior strongly suggests d-wave pairing, with line nodes. 

Measurements of Tf^ are straightforward in the CeCoIns because the In is spin 
/ = I and has a large quadrupolar frequency. The data of several groups show consistent 
behavior, consistent with line nodes, as seen in Fig. [HI Measurements of T^^ in 
field have not been performed, but may be important since recent small angle neutron 
scattering data suggest the presence of antiferromagnetic vortex cores in this material, 
which would contribute extra relaxation channels [87]. This antiferromagnetism is 
probably related to the proximity to a quantum critical point. 

6.2.4- PuCoGa^ and the bridge to high Tc superconductivity Until the year 2002, 
the class of unconventional superconductors seemed to consist of three separate 
subclasses: the high-Tc cuprates, the organic superconductors, and the heavy fermion 
superconductors. The cuprates, with transition temperatures around 100-150 K, and 
the organic superconductors, with transition temperatures on the order of 5-lOK, have 
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Figure 14. T[' versus temperature in CeCoIus measured by NQR at zero applied 
field, reproduced from [66]. The solid black line is a fit to T^, as expected for line 
nodes in the superconducting gap. 




Figure 15. The Knight shift of the Co in PuCoGas, showing the suppression of the 
spin susceptibility below Tc, indicative of spin singlet pairing. 



NMR in heavy fermion superconductors 



35 



structures based on two dimensional doped Mott insulators. The heavy fermions 
have transition temperatures of 1-2 K, are good metallic systems with either Ce or 
U atoms in the lattice. Since the physics of the basis materials seem to be quite 
different, it is natural to suspect the nature and origin of the superconductivity may be 
different as well. However, these different sets of materials have similar phase diagrams, 
with superconductivity emerging near to antiferromagnetism, and in most cases the 
superconductivity appears to be unconventional. In 2002, PuCoGas was found to be 
superconducting at 18.5 K. This result was surprising because the transition temperature 
is so high and the structure is identical to that of CeCoIus. By going from the 4f to 
the 5f series, the transition temperature increased by an order of magnitude. One of 
the leading questions was whether the superconductivity in PuCoGas is conventional or 
unconventional. NMR played a crucial role, as one of the only techniques available to 
deal with the radioactive plutonium compound. Knight shift measurements of the Co 
and Ga clearly revealed a spin singlet superconductor, and T^^ measurements showed a 
nodal gap. The fact that PuCoGas was a heavy electron unconventional superconductor 
suggested then, that the gap between the heavy fermions and the high-T^ cuprates may 
be artificial, and that unconventional superconductivity can arise over a continuous 
range of temperatures, possibly controlled by the scale of antiferromagnetic exchange in 
the material [26l[63]. Another Pu-based superconductor, PuRhGas with Tc ~ 9 K, was 
discovered shortly after PuCoGas, and NMR measurements of K and T^^ also reveal an 
nodal, spin-singlet superconductor [88]. Quite recently in 2007, another actinide-based 
superconductor was discovered, NpPd5Al2 (Tc = 4.9 K) [H9l|90]. NMR measurements 
reveal line nodes and spin-singlet behavior in this system, also suggestive of d-wave 
superconductivity [91]. Clearly, the actinide series holds many surprises to come, and 
NMR will play a critical role in determining the properties of any superconducting 
materials. 

There are other heavy fermion superconductors, in particular UBeis, UPd2Al3, 
UNi2Al3, where NMR measurements in the superconducting state have shown power- 
law behavior of the spin lattice relaxation rate in the superconducting state, but the 
Knight shift measurements and overall interpretation have remained inconclusive. We 
refer the reader to [92] for more details. Other materials, such as CeRhlns and Celns 
become superconducting under pressure. In these cases, T^^ measured by NQR reveals 
behavior, but depending on the pressure there is evidence for a T— linear term at 
lower temperatures [931 [Ml [95] . This crossover behavior at low temperature is poorly 
understood, but may be related to the coexistence of antiferromagnetism [96] . 

6.3. The field-induced phase in CeCoIn^ 

Some of the most striking and unusual behavior exhibited by CeCoIns emerges at 
ambient pressure in the superconducting state. In 2002, heat capacity measurements 
close to Hc2 revealed an unexpected first order transition between the normal state 
and the superconducting state, and further measurements indicated the presence of a 
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Figure 16. The superconducting phase diagram of CeCoIns (structure shown as 
inset). The solid black hues represent second order phase transitions and the solid 
blue line represents a first order transition. The A phase is a conventional Abrikosov 
vortex lattice, and the B phase is the field-induced phase that may be related to an 
FFLO phase. The phase boundary between the A and B phases is referred to in the 
text as either T*{H) or H*{T). 

second order phase transition within the superconducting state a temperature T*{H) 
[971 EH El]. This phase was initially identified as a possible Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) phase, which was postulated to exist in superconductors with 
exceptionally large spin susceptibility [1001 IIOH I102j . Magnetic fields can suppress 
superconductivity by two mechanisms: orbital limiting or Pauli limiting. In an orbital 
limited superconductor, the highest magnetic field that the system can support is limited 
by the critical current density, Jc- The London equations show that superconductors 
respond to external fields by creating supercurrents to screen the internal field, but if the 
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supercurrent exceeds Jc then the system can no longer support superconductivity [IS] . 
In a Pauh-hmited superconductor, the maximum magnetic field is determined by the 
spin susceptibility of the quasiparticles. For sufficiently large fields, the Zeeman energy 
of the quasiparticles can exceed the superconducting condensation energy, in which case 
the material can no longer support superconductivity. Typically the Pauli susceptibility 
of superconducting materials is small enough that the limiting field is determined by 
the orbital limit. However, since heavy fermions typically have large susceptibilities due 
to the enhanced effective mass, the Pauli limit may be reached in principle. 

In Pauli limited superconductors, the transition between the normal and 
superconducting phases can become first-order, and in some cases a new phase can 
emerge, in which the superfiuid density becomes modulated over large length scales. 
In the FFLO phase, the order parameter \E' ~ \E'oexp(— iq ■ r) or \E' ~ \E'oCOs(q ■ r), 
where q = kp^^i — ^f,i, where ^f,'[,1 Fermi wavevectors of the spin-up and spin- 

down Fermi surfaces. The wavelength, A = 27r/g, is therefore several times the unit 
cell length, a. The minimum possible distance must be Amin = ^, the superconducting 
coherence length, since the condensate cannot respond on a length scale shorter than 
^. In this phase, the wavefunction vanishes at nodes in real space, and the unpaired 
up-spin quasiparticles reside in these nodes. This phase was predicted to exist nearly 
forty years ago, but experimental evidence for its existence has been sparse. Likely 
candidates include the organic superconductor (TMTSF)2C104 and CeCoIns |103[ [99] . 
The strongest evidence for the existence of this phase in CeCoIus stems from the 
similarity of the phase diagram, shown in Fig. [161 to calculations. All of the theoretical 
calculations, however, are based on either a non-interacting or weakly interacting Fermi 
liquid [1041 1105] . The CeCoIns, however, clearly has more physics at play than a simple 
Fermi liquid, as there are strong f-electron exchange interactions (RKKY) and Kondo 
screening of the moments, leading to non-Fermi liquid behavior in the normal state 
[106[ I107j . Therefore, it is not obvious a priori that this new phase can be described 
by the simple theories of FFLO developed for Fermi liquids. Indeed, NMR evidence 
suggests that the situation is more complex. We therefore refer to this phase simply as 
the B phase, in order to distinguish it from the classical FFLO phase described by these 
theories. 

The first NMR experiments in the B phase were reported by ^108| . The spectra 
of the In(l) in this report showed dramatic changes as the temperature was lowered 
through T*. Below T*, the spectra revealed two components, one peak at slightly 
reduced frequency, and a second smaller peak at the resonance frequency of the normal 
state (see Fig [17]). These data were interpreted as evidence for coexistence of both 
superconducting and normal regions, consistent with the FFLO picture of normal planes 
sandwiched between superconducting layers. One puzzling aspect of this data, however, 
is that the spectra seemed to reveal a continuous evolution rather than a first order 
transition. Secondly, the change AK = K{Tc) — K{T = 0) of the Knight shift of 
the lower (superconducting) peak was much smaller than AK reported previously at 
lower fields [TOl [H] . Subsequent measurements by other groups at high fields disagreed 
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Figure 17. The Knight shift of the In(l) in CeCoIns for H\\ab as measured at various 
field and reported by Kakuyanagi et al. [108j . Young et aL [SI], and Mitrovic et al. 
[52] . Dotted hues are guides to the eye for discontinuous transitions. 



with this initial report, both quantitatively and in interpretation of the data 
The data sets are compared in Fig. [T7J Mitrovic et al. argued that the spectra not 
only showed much larger AK, but also that the intensity of the signal is suppressed by 
roughly an order of magnitude below T^. The original Kakuyanagi et al. data showed 
no decrease in intensity below Tc, although in practice the signal is suppressed because 
the rf penetration is reduced. Young et al. reported discontinuous jumps in the Knight 
shift, with significantly larger AK values. As seen in Fig. [TTJ the Young and Mitrovic 
data agree well with one another, but not with that of Kakuyanagi. Recently, Mitrovic 
et al. published a comment on the original Kakuyanagi paper, in which they argue 
that the Kakuyanagi data are incorrect and are the extrinsic effects of rf heating of the 
sample [109] . Mitrovic et al. were able to reproduce the two peaks observed in the 
Kakuyanagi data by systematically increasing the rf power. Since the transition is first 
order, it is reasonable to expect that there may be some coexistence of both normal 
and superconducting regions near Tc, and for excessive rf heating (Joule-heating by 
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Figure 18. (Left panel) NMR spectra in tlie normal state {T — 0.475 K) and in the 
B phase (T = 0.050 K) of CeCoIng at 11.1 T, showing the spectra of the Co (/ = 7/2), 
several of the In(l) (/ = 9/2) transitions, and one of the In(2) transitions. Clearly, 
several differences emerge in the B phase, most notably the broad double-peak form of 
the In(2), which is characteristic of an incommensurate magnetic order. In contrast, 
the In(l) and Co show little change. (Right panel) Similar spectra, but at a higher 
field of 11.485 T, in the normal state (T = 0.36 K) and in the B phase (T = 0.05 K). 

induced currents in the metallic sample) this coexistence may persist over much larger 
temperature ranges. 

Young et al. also report data on the Co and In(2) sites over a series of temperatures 
crossing from the normal state, through the mixed phase, and into the field-induced 
phase. They found that the Co, the In(l) and the In(2) resonances behaved very 
differently from one another in both the mixed phase and the B phase. As seen in 
Fig. [181 the In(2) reveals a broad two-peak spectrum consistent with antiferromagnetic 
order. The resonances of the other two sites, on the other hand, remain relatively sharp. 
Young et al. argued that since the response of the three sites in the unit cell could not be 
explained by the conventional picture of the FFLO phase, in which the inhomogeneity is 
a long-range phenomenon with A ^ a, that a more complex scenario must be present. 
In particular, the B phase exhibits long range antiferromagnetism. The spectrum of 
the In(2) in this phase is similar to that of the In(2) in the antiferromagnetic state 
of CeRhlus, and CeCo(Ini_2;Cd2:)5 |110j . Assuming isotropic hyperfine couplings, the 
NMR data can be explained by an incommensurate structure with Q = (| -|- 5, |, |) 
and moments along the applied field direction, which gives rise to a finite hyperfine 
field at the In(2) sites, but not at the In(l) or Co sites, as observed in the data. This 
wavevector represents only one possible magnetic structure that can explain the set of 
hyperfine fields observed in the B phase and is not unique, but is the most likely one if 
the hyperfine couplings to the In(2) are purely isotropic. 

Recently elastic neutron scattering experiments were reported in the B phase of 
CeCoIns, which clearly reveal an incommensurate magnetic order and thus confirming 
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the work of Young and coworkers In these experiments, the field was ahgned 

along [1-10], and the Ce moments were observed to lie along [001] with a magnitude 
of 0.15/xb, and wavevector Q = (| + 5, | + 5, |) with 6 = 0.06. In this case, the 
incommensurate modulation lies in the [110] direction, perpendicular to both the applied 
field and the ordered moments. These observations contrast somewhat with the original 
NMR measurements with field along [100]. Although it is possible that the ordering 
wavevector and the direction of the moments may change when the applied field is 
rotated by 45°, the discrepancy between the two measurements probably arises from 
insufficient information about the hyperfine tensor of the In(2). In the NMR report, the 
tensor was assumed to be isotropic, but if this constraint is relaxed, then other magnetic 
structures consistent with the spectra of the In(2) emerge as possible candidates. In 
particular, if the tensor has dipolar symmetry, then ordered moments along [001] can 
give rise to the observed spectra. 

One may then ask whether the B phase is an FFLO phase or an antiferromagnetic 
one. In fact, identifying this phase as an FFLO phase is nearly impossible with NMR, 
since one would have to detect a change of phase of the order parameter along the spatial 
modulation direction. However, bulk measurements of the phase diagram are consistent 
with theories for an FFLO phase in a Pauli limited superconductor [105[ 11041 1102[ I112j . 
Furthermore, NMR measurements clearly show no sign of any magnetism for fields 
H > Hc2- If the B phase were a secondary phase that either competes with, or 
coexists with the superconductivity, then it should persist above Hc2 in the normal state. 
Indeed, in CeRhlus under pressure, there is a field-induced phase that coexists with 
superconductivity below Hc2 and exists on its own above (see Fig.) [113]. In contrast, 
the B phase of CeCoIns exists only within the superconducting phase. This might be 
explained by a subtle change in the delicate balance between the RKKY interaction and 
Kondo screening of the Ce moments. Indeed, the onset of a superconducting gap can 
modify both these interactions, and may have a stronger effect on the Kondo screening 
|114[ [Mj. Therefore, for H ~ Hc2, where the superconducting condensation energy, the 
RKKY interaction and the Kondo interaction are all of the same magnitude, it might be 
possible that the Kondo interaction is suppressed in the presence of the superconducting 
gap, and the balance is tipped in favor of the magnetic order. Above Hc2, when the gap 
vanishes, the balance is reversed, and magnetism would disappear. 

An alternative explanation is that the driving element for the B phase is truly 
an instability towards the FFLO state, and the the antiferromagnetism grows in as 
a secondary parameter that may nucleate in the normal planes of the FFLO phase. 
This scenario offers an explanation for the pressure dependence of the superconducting 
phase diagram, which shows that with increasing pressure the volume of phase space 
occupied by the B phase increases |115j . Pressure tunes the CeCoIus away from 
the antiferromagnetic quantum critical point, and suppresses the antiferromagnetic 
fluctuations [65l 166] . If the B phase were driven by an antiferromagnetic instability, 
one would expect the opposite trend with pressure. In fact, there is evidence even 
within the mixed phase for antiferromagnetism within the vortex cores of the Abrikosov 
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vortex lattice [MlEZl- Small angle neutron scattering (SANS) data for the field along 
the c-direction indicate that the form factor for the vortex lattice increases with field, 
rather than decrease as expected for a conventional vortex lattice [37] . Furthermore, 
the NMR linewidths of the Co, In(l) and In(2) in the mixed phase are significantly 
larger than expected by Eq. (l28l) . In fact, the second moment depends on the 
particular nucleus, which cannot be explained if the broadening is purely due to orbital 
magnetisation. Rather, there must be a spin magnetization contribution as well, which 
will give rise to different hyperfine field distributions at the different nuclei. This 
inhomogeneous spin distribution is probably due to local moments in the vortex cores. 
Bulk magnetization measurements also indicate a strong paramagnetic contribution 
in the mixed phase, which is further evidence for magnetism in the cores |116j . A 
secondary antiferromagnetic order parameter that nucleates in the vortex cores where 
the superconducting order parameter goes to zero is a natural expectation of a Ginsburg- 
Landau description of competing order parameters, and is also observed in the high-Tc 
cuprates |117l 1118] . If the superconducting order vanishes at nodal planes, rather than 
simple vortex lines, then it may be easier for a competing order parameter to establish 
long range order. 

Importantly, the neutron scattering data of Kenzelmann et al. reveal that the 
ordering wavevector Q is independent of applied field in the B phase. This result 
is inconsistent with an FFLO scenario, where the instability is driven by the Pauli 
susceptibility. On the other hand, using a Landau phenomenological theory for the 
coupling between the field and the superconducting order parameter, Kenzelmann et 
al. argue that the most likely couphng term leads to a superconducting gap function 
that carries finite momentum. It appears, then, that this phase consists of a pairing of 
electrons with different wavevectors, as in an FFLO phase, but the driving interaction 
may be an antiferromagnetic instability instead of the large Zeeman term that leads to 
a conventional FFLO picture. Clearly, this new phase offers a wealth of new physics, 
and NMR is likely to play a central role in its study. 

6.4- Probing the coexistence of antiferromagnetism and superconductivity 

NMR has played a central role in studies of coexisting antiferromagnetism and 
superconductivity in a number of heavy fermion compounds. These two order 
parameters may arise in different channels, for example if there are more than one 
type of magnetic atom in the unit cell |119j . This may be relevant for several U-based 
heavy fermions, including UPd2Al3, UPts, URu2Si2, and UNi2Al3, since the U atom 
in these cases is typically has more than one 5f electron. Here we focus on studies 
of coexisting orders for Ce compounds, in which there is only one channel available 
for the 4f^ electronic degree of freedom from the Ce^+ ions. The same degree of 
freedom is responsible for the local moment magnetism, the heavy fermion state, and 
the unconventional superconductivity. This is clearly seen in NMR studies of CeRhlns 
under pressure |9ll [931 1120] . The temperature-pressure phase diagram of CeRhlus is 
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shown in Fig. (1191) . For pressures between ~ 1.5 — 1.8 GPa, and T < the system 
shows bulk superconductivity microscopically coexisting with antiferromagnetism, as 
seen in the NQR spectra of the In(l). The spectra clearly show that all of the 
In(l) nuclei in the same experience a static hyperfine field (see §2.5). If the material 
consisted of an inhomogeneous coexistence of regions that are superconducting but 
not magnetic, and regions that are antiferromagnetic but not superconducting, then 
the spectra would show two sets of resonances: one that is split by the internal field 
corresponding to the antiferromagnetic regions, and one that is not split, corresponding 
to the superconducting regions. In other words, the spectra would be a combination of 
both the normal state spectrum and the ordered state spectrum as shown in Fig. (ITQll . 
Clearly, this is not the case. Furthermore, T^^ measurements clearly show anomalies at 
both Tn and Tc. Although this observation is not as strong an argument for coexistence 
as the spectral evidence, it does support the interpretation since presumably the same 
set of nuclei experience both phase transitions. It should be pointed out, however, 
that T]~^ is measured by fitting the relaxation of the bulk magnetization signal, which 
can in principle be composed of two signals: M{t) = Mi{t) + M2{t), where each set 
experiences only one of the two phase transitions. In this case, the magnetization 
recovery will show multiple relaxation time scales. Therefore, it is crucial to verify 
whether the magnetization recovery shows a single or multiple relaxation time scales 

m- 

In contrast to the CeRhlus under pressure, NQR experiments in Celus reveal 
antiferromagnetism and superconductivity do not coexist microscopically, but are phase 
separated ^95j. Indeed, the In NQR spectrum consists of a superposition of the static 
broad line associated with the antiferromagentic order, and a narrow line associated 
with paramagnetism (the superconducting regions). The relative volume fractions of 
these two regions vary as a function of pressure, with antiferromagnetism giving way 
to superconductivity, and the maximum in Tc occurs when the two volume fractions 
are the same. In both the CeRhlus and the Celua case, the entropy associated with 
the two phase transitions arises entirely from single Ce 4f electron. It is not clear why 
this degree of freedom can exhibit multiple orders simultaneously in the case of the 
CeRhlns, and not in the Celna, but may be related to the presence of a first order phase 
transition in the latter case [95]. As a third example, in CeCu2Si2, superconductivity and 
antiferromagnetism coexist microscopically in some portions of the phase diagram, and 
in other regions the antiferromagnetism abruptly disappears when superconductivity 
sets in [93 |ll23|,ll24j . In the latter case, a possible scenario is that the antiferromagnetism 
is a spin density wave (SDW) order, and both order parameters compete for regions of 
the Fermi surface. In contrast, both the fluctuations and the phase diagram in CeRhlus 
suggest that a Fermi surface competition is not the case, but rather a dramatic Fermi 
surface reconstruction occurs in which the Ce 4f electrons change from localized to 
itinerant character [122] . 

It is interesting to note that in states with coexisting antiferromagnetism and 
superconductivity T^^ typically exhibits a crossover in behavior in the superconducting 
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Figure 19. The temperature-pressure phase diagram of CeRhlns at zero field and 
at 3.3 T (reproduced from [113]). The spectra of the In(l) (±1/2 ^ ±3/2) NQR 
transitions (from [120j ) are shown at the indicated points in phase space, at = 0. 
In the range between 1.5 - 1.8 GPa, antiferromagnetism and superconductivity coexist 
microscopically. In field, the region of coexistence extends up to higher pressure up to 
2.1 GPa (not shown) where deHaas-van Alphen experiments indicate a change of the 
Fermi surface |122j . 



state. This is seen in CeRhlns under pressure |9l], in CeRhi_a;Ir^In5 [96], and in 
CeCo(Ini_^Cd^)5 [ml [125]. In CeRhlug and CeRhi.^Ir^Ing where Tj^ > T„ Tf ^ 
below Tc initially drops as T^, then crosses over to T at some temperature T < T^. 
This might be understood as some portion of Fermi surface that remains ungapped. In 
the CeCo(Ini_a;Cd2:)5, the low temperature behavior behaves as T^/*^, but it is not clear 
why this should be so. However, if the antiferromagnetism and superconductivity coexist 
microscopically, then it is clear that there are ordered local moments that create a field 
at the nuclear sites. Goldstone fluctuations of the antiferromagnetic order can give rise 
to fluctuating hyperfine fields at the nuclear sites that potentially can relax the nuclei. 
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At present, it is not clear if the mechanism of the spin lattice relaxation rate in this state 
arises from these magnons or from itinerant quasiparticles from gapless regions of the 
Fermi surface. Evidence for the former scenario can be found in comparisons of T^^ at 
the In(2) versus In(l) site in CeRhi.j-Ira-Ins, which show an onset temperature for the 
T-linear behavior higher for the In(2) than for the In(l) [96]. Since the static hyperfine 
field at the In(2) is higher, it is reasonable that the T['J;^^g^^{ln{2))> T{;^„g„„„(In(l)), 
and therefore becomes the dominant relaxation mechanism at a higher temperature. If 
the relaxation mechanism were purely from excited itinerant quasiparticles, then the 
temperature dependence of the In(l) and In(2) would be identical. 

6.5. Impurities in Heavy Electrons 

6.5.1. Radioactive decay in PuCoGa^ Superconductivity responds to impurities in 
a predictable fashion that depends on the nature of the superconductivity |126i 1^ . 
Since NMR is sensitive to the density of states and is a microscopic probe, it 
is an ideal technique for investigating the nature of localized impurity states in 
superconductors and correlated electron systems |127] . In a d-wave superconductor, 
pair breaking impurities locally suppress the superconducting order parameter and give 
rise to localized electronic states with energies below the superconducting gap. With 
increasing impurity concentration, these localized states can overlap and develop into an 
impurity band that can contribute to the spin lattice relaxation rate. In many d-wave 
superconductors, this impurity band short-circuits the behavior and for sufficiently 
low temperature, crosses over to T-linear behavior |128j . A clear example of this 
behavior is seen in PuCoGas, where the radioactive decay of the Pu strongly affects the 
superconductivity. As ^^^Pu— 235^ _j_ 4-|^g j^^^g ^ YiqX{ life of 24,110 years, significant 
lattice damage can build up in this material over time as the recoil of the heavy U 
nucleus, with a kinetic energy of 86keV, displaces several thousand atoms and the 
energetic alpha particle deposits its energy in the lattice |129] . For sufficiently high 
temperatures, the lattice damage and Frenkel pairs can anneal out |130] . Over time, 
however, pair breaking impurities suppress the superconductivity and is reduced at 
the rate of dTc/dt ^ 0.2 K/month. Microscopically, these pair breaking impurities give 
rise to an impurity band near Ep and prevent the density of states from reaching zero 
at Ep as it would in a pristine superconductor. This increase in iV(0) is reflected as an 
increase in T^^ and gives rise to a linear term that becomes important at sufficiently 
low temperature, and a finite as T — >^ 0, rather than — >^ as T — 0. This effect is 
clearly observed in Knight shift and measurements in PuCoGas; indeed the size of 
these effects correlates well with the age of the sample and the decrease of with time 

m- 

6.5.2. Magnetic droplets in quantum critical systems Clearly, NMR can shed important 
light on the nature of the superconductivity by probing how the local density of states 
is modified by the presence of impurities. Another clear example in which NMR can 
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elucidate the nature of the physics associated with impurities is in doped quantum 
critical systems |13H I132[ 11331 I134j . In a system that is close to a quantum phase 
transition, impurities may induce local perturbations that nucleate droplets of the 
ordered phase. NMR in such a system may be sensitive to this local order. In 
Cd doped CeCoIns, the Cd nucleates antiferromagnetic order with radius C,af- For 
sufficiently high doping, these droplets interact with one another giving rise to long range 
antiferromagnetism |135] . The Cd doping also suppresses superconductivity, and above 
a critical concentration of dopants the system is purely antiferromagnetic. The dynamics 
of this change in behavior is clearly visible in the T^^ of the In(l). Furthermore, since 
the Cd dopants nucleate antiferromagnetism, the local spin magnetisation develops a 
staggered response, similar to the Friedel oscillations induced by impurities in a Fermi 
liquid [127] . The Co in this system is particularly sensitive to this staggered response, 
since it is located in a symmetric site between nearest neighbor Ce atoms. The Co 
experiences a transferred hyperfine interaction to these two Ce spins, so that the form 
factor at the Co site vanishes for antiferromagnetic correlations. The hyperfine field 
vanishes, however, only if the two Ce spins are equal in magnitude. Around an impurity, 
the magnetization will decay within the correlation length (the details of how this order 
varies spatially depend critically on the nature of the quantum critical system) [133j . 
Therefore, Co atoms sitting within such a droplet will experience finite hyperfine fields. 
This distribution of fields is manifest as a broad linewidth that varies with doping and 
temperature (see inset Fig. [201) . For non-interacting droplets, one can show that the 
linewidth. A, is given by + A^^(T) where Acd(7') ~ McdiT)^AF{T). Here Aq is 
the intrinsic Co linewidth (Eq. [5]), and Mcd(r) is the magnetization of the droplets. 
Fig. [20] shows Acd(T) / Mcd(T) as a function of temperature, which clearly shows that 
^AF increases below ~ 10 K. A complete analysis is not available in this case, however, 
since the Cd concentration [x =1.5%) is high enough that the droplets interact. Similar 
effects were observed in the Ni and Zn doped cuprates |136l 11371 11381 1139] . This new 
method promises to shed new light on the dynamics of the antiferromagnetic correlations 
that develop in these and other heavy electron materials. 

6.6. Summary 

At the time of this writing, a new class of iron oxypnictide superconductors has just 
been discovered, and NMR is playing a central role in experimental studies of the 
superconductivity and antiferromagnetism in these compounds. Clearly, NMR has 
the potential and versatility to shed important light on the physics of these systems, 
as well as other classes of materials that may be discovered in the future. There 
remain still several open questions in the study of the heavy fermions, though, that 
can be addressed with NMR. In particular, the interplay between antiferromagnetism 
and superconductivity in different classes of materials and how these order parameters 
emerge for different classes of quantum critical behavior remains a fascinating area of 
research. Secondly, the behavior of non-centrosymmetric superconductors is an active 
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Figure 20. A plot of Acd(r)/Mcd(r) - ^ A^T) - A''J (xcdiT) - xo(T)) ~ ^sf 
versus temperature in CeCo(Ino.85Cdo.i5)5. Tlie Inset shows the FWHM, A(T), of the 
Co resonance in the same material. Clearly, the antiferromagnctic correlations begin 
to grow in below ~10 K. 



area of research for NMR [140j . And finally, the technique of purposefully adding 
impurities to these strongly correlated systems as a means to investigate the local 
spectroscopy of the impurity states with NMR holds tremendous promise to shed light 
on the nature of the order parameters that can emerge. 
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